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Dear Sir: 



FIRST PRELIMINARY AMENDMENT 

Please enter the following preliminary amendment to the application: 



Replacement Specification 

Please delete the specification currently on file and replace with new specification, pages 
1 through 21, as submitted in Attachment A. 

In the Claims 

Please delete claims 1 through 28 currently on file and replace with new claims 1 through 
28, as submitted in Attachment A. 

In the Abstract 

Please delete the page containing the Abstract currently on file and replace with the new 
page containing the Abstract, as submitted in Attachment A. 

In the Drawings 

Please replace Figures la, lb, Ic, 2a, 2b, 3a, 3b, 4a, 4b, 5, 6a and 6b currently on file 
with new Figures la, lb, Ic, 2a, 2b, 3a, 3b, 4a, 4b, 5, 6a and 6b, as submitted in Attachment A. 

Attachment A is a replacement specification, claims, abstract and drawings. 
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Attachment B is a clean copy of the amendment. 



Attachment C is a marked up copy of the amendment. 



REMARKS 



To facilitate the examination process, a replacement specification, claims, abstract and 
drawings is provided as Attachment A. No new matter has been added. This document is 
identical to the international publication number WO 01/53982 Al. 

We request a statement be added at the beginning of the specification under 35 U.S.C. 
1 19(e) and in accordance with 37 CFR 1 .78. A clean copy of the new paragraph and a marked up 
copy of the new paragraph is attached in Attachment B and Attachment C. 

Applicant looks forward to receiving the first Official Action from the Examiner in due 

course. 

Respectfully submitted. 




Gordon Freedman 
Reg. No. 41,553 



Freedman & Associates 
117 Centrepointe Drive 
Suite 350 
Nepean, Ontario 
Canada K2G 5X3 



Tel (613) 274-7272 
Fax (613)274-7414 
Email: gordon@ipatent4u.com 
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ATTACHMENT B 

Clean New Paragraph 

After the Title 

Page I, please insert new paragraph 000 directly after the title: 

[000] This application is a national phase entry of International Application No. 

PCT/CAO 1/0003 8 filed January 17, 2001, which claims priority firom United States Provisional 

Patent Application No. 60/176,298 filed January 18, 2000, now abandoned. 
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ATTACHMENT C 

Marked Up New Paragraph 

After the Title 

Page 1, please insert new paragraph 000 directly after the title: 

[OOP] This application is a national phase entry of International Application No. 
PCT/CAOl/00038 filed January 17. 2001. which claims priority from United States Provisional 
Patent Application No. 60/176.298 filed January 18. 2000. now abandoned. 
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PARALLEL CORRELATOR ARCHITECTXTRE 



Field of the Invention 

The invention relates generally to spectral analysis of wide-band sampled signals, 
5 in particular to a method and system for real time digital spectral analysis of wide-band 
sampled signals using parallel processing techniques. 

Background of the Invention 

SpKJtral analysis is a well-icnown analytical tool, which has been used for decades 
in science and industry to characterize materials, processes and systems based on spectral 

10 information. For example, cross-correlation of spectral data provides amplitude and phase 
information between two signals, which have some common component generated by 
some common artificial or natural process. In radio astronomy, cross-correlation allows 
very weak signals from distant astronomical sources to be extracted from strong 
independent system noise. Amplitude and phase from many different antennae are then 

15 used to build an image of the radio source. In network analysis an artificial signal is 

generated, transmitted through an object under test, and is then cross-correlated with the 
original signal. This yields the amplitude and phase response as a function of frequency 
of the tested object. In aperture synthesis radar an object is irradiated by a high power 
microwave source. Cross-spectral analysis of the microwave reflections from the object's 

20 surface is used to produce an image of the object. Auto-correlation of a signal is used in 
spectrum analysis to find periodic signal components in signals disturbed by random 
noise. Another field for the application of spectral analysis is in phased array systems 
such as phased array radar, remote sensing phased array receivers and wide-band phased 
antennae. 

25 In most of these applications it is desired to use real time digital signal processing 

to determine the spectrum of a sampled analog signal. Known signal processing systems 
use mixed analog and digital systems. Furthermore, in current phased array systems 
analog methods for sub-sample delay interpolation are used. However, analog systems 
are inherently unstable due to time and temperature variability of analog filters. 
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Therefore, real time digital processing systems using time-domain de- 
multiplexing have been recently developed. Unfortunately, these systems require more 
than an order of magnitude more digital hardware for large array systems than a direct 
method, as noted helow. Furthermore, these systems are often limited by clock rates of 
5 the digital equipment that are much lower than the bandwidths of signals to be processed. 
It is desirable to use digital techniques in real tune to determine the spectrum of a 
sampled analog signal whose bandwidth exceeds the capability of realized digital systems 
to obtain the spectrum directly. A dh-ect method according to the prior art would require 
digital processing to occur at clock rates of twice the bandwidth of the signal. 

10 Recent developments in computer technology overcome the limitations of 

processor clock rates by processing one task on a plurality of processors in parallel in 
order to process a large amount of data processed in a fraction of the time needed by a 
one-processor system. 

It is, therefore, an object of the invention to provide a method and system based 
15 on parallel processing techniques to determine a wide-band spectrum efficiently with 

arbitrarily high spectral resolution using a processor clock rate that is an arbitrary fraction 
of the wide-band clock rate. 

It is another object of the invention to provide a precision sub-sample digital delay 
interpolation for phased array and cross-correlation systems. 



Summary of the Invention 

In accordance with the present invention there is provided a method for real-time 
digital spectral analysis of wide-band signals comprising the steps of: 
25 receiving a wide-band signal; 

shifting the center frequency of the wide-band signal by a small fraction e of its 
bandwidth; 

sampling and digitizing the wide-band signal; 
processing the digitized wide-band signal using a digital filter; and, 
30 decimating the digitally filtered wide-band signal. 



20 
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In accordance with the present invention there is further provided a method for 
real-time digital spectral analysis of wide-band signals comprising the steps of: 
receiving a wide-band signal; 

shifting the center frequency of the wide-band signal by a small fraction f of its 
5 bandwidth; 

sampling and digitizing the vdde-band signal; 

de-multiplexing the digitized wide-band signal into N parallel sample streams; 
processing the 7v^ parallel sample streams in parallel using N digital FIR filters; and, 
determining 2'' •N;k = 0,1,... sub-band signals by decimating the sample stream from 
10 each FIR filter by a factor of 2* • 7/; A: = 0,1,.„, wherein only every 2* • N"';k = 0,1,... 
sample is retained and the others are discarded. 

In accordance with an aspect of the present invention there is provided a method 
for cross-correlating de-rotated sub-band signals sub-band by sub-band, the method 
15 comprising the steps of: 

receiving 2* • N;k = 0,1,... pairs of first and second de-rotated sub-band signals 
at2^ ■ N;k = 0,1,... cross-correlators, wherein each pair is received at a different cross- 
correlator of the!'' -Nik - 0,1,... cross-correlators; 

delaying one of the first and second de-rotated sub-band signals with respect to the other 
20 in a series of delay intervals at each of the 2* - JV; A: = 0,1, ... cross-correlators; 

forming the product of the first and the second de-rotated sub-band signals at each of the 
delay intervals at each of the 2* ■ N;k = 0,1,... cross-correlators; 
producing a sub-band cross-correlation result at each of the 2* -Nik- 0,1,.-. cross- 
correlators by summing the products over a period of time; 
25 transforming each sub-band cross-correlation resuU at each of the 2* ■ Nik = 0,1,... cross- 
correlators by means of a Fourier Transform into a cross-spectrum result; and, 
correcting each cross-power spectral point of each sub-band cross-spectrum result with a 
sub-band scaling term, a gain differential compensation term, a bandshape correction 
term and a wide-band power gam term. 
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In accordance with the aspect of the present invention there is further provided a 
method for real-time digital j pectral analysis of wide-band signals comprising the steps 
of: 

receiving a first and a secona wide-band signal; 
5 shifting the center frequency of each of the first and the second wide-band signal by a 

small fraction s, and , respectively, of its bandwidth; 

sampling and digitizing the first and the second wide-band signal; 

de-multiplexing each of tiie digitized first and second wide-band signals into first AT 

parallel sample streams and second N parallel sample streams; 
10 processing each of the first and the second N parallel sample streams in parallel using 

2- AT digital FIR filters; 

determining first and second sub-band signals by decimating the sample stream from 
each FIR filter by a factor of , wherein only every N"" sample is retamed and the 
others are discarded; 

15 re-quantizating the N sub-band signals by re-scaling and truncating in order to reduce 
downstream processing load; 

phase rotating each of the first and second N sub-band signals by phase £j and > 
respectively, using a digital phase rotator producing first and second JV de-rotated sub- 
band signals; 

20 receiving pairs of the first and second N de-rotated sub-band signals atiV cross- 
correlators, wherein each pair is received at a different cross-correlator of theiV^ cross- 
correlators; 

delaying one of the first and second de-rotated sub-band signals with respect to the other 
in a series of delay intervals at each of theiV cross-correlators; 
25 forming the product of the fu^ and the second de-rotated sub-band signals at each of the 
delay intervals at each of the 7^ cross-correlators; 

producing a sub-band cross-correlation result at each of the N cross-correlators by 
sunmiing the products over a period of time; 

transformmg each sub-band cross-correlation result at each of the N cross-correlators by 
30 means of a Fourier Transform into a cross-spectrum result; 

correcting each cross-power spectral point of each sub-band cross-spectrum result with a 
sub-band scaling term, a gain differential compensation term, a bandshape correction 
term and a wide-band power gain term; and. 
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concatenating flie sub-band cross-spectra to obtain a wide-band spectnim. 

In accordance with ariother aspect of the present invention there is provided a 
method for real-time digital spectral analysis of wide-band signals of a phased array 
5 system comprising the steps of: 

receiving M wide-band signals from the phased array system; 

shiftmg the center frequency of each of the M wide-band signals by a small fraction 

,...,s„ , respectively, of its bandwidth; 
sampling and digitizing the M wide-band signals; 
10 de-multiplexing each of the digitized M wide-band signals into M-N parallel sample 
streams; 

processing each of the M-N parallel sample streams in parallel using M-N digital FIR 
filters; 

determining M- N sub-band signals by decimating the sample stream from each FIR 
15 filter by a factor of N , wherein only every N"" sample is retained and the others are 
discarded; 

re-quantizing the M- iV sub-band signals by re-scaling and truncating; 
complex mixing each of the M-N re-quantized sub-band signals; 
90** phase shifting one of two components of each of the complex mixed M- N sub-band 
20 signals; and, 

forming multiple beams by adding same sub-band components of the M-N sub-bands, 
the same sub-band components being provided by same digital FIR filters. 

In accordance with the present invention there is further provided a system for 
25 real-time digital spectral analysis of wide-band signals comprising: 
a port for receiving a wide-band signal; 

means for shifting the center frequency of the wide-band signal by a small fraction s of 
its bandwidth; 

an A/D converter for sampling and digitizing the wide-band signal; 
30 a de-multiplexer for de-multiplexing the digitized wide-band signal into N parallel 
sample streams; and. 
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N processors, each processor processing one of the N parallel sample streams by 
digitally FIR filtering and decimating the sample stream in order to determine a sub-band 
signal. 

5 Brief Description of Figures 

Exemplary embodiments of the invention will now be described in conjunction 
with the following drawings, in which: 

10 Figure 1 a is a simplified block diagram of a method according to the invention for 
spectral analysis of a wide band signal; 

Figure lb is a simplified flow diagram of a method according to the invention for spectral 
analysis of a wide band signal; 

15 

Figure Ic is a simplified block diagram of a system according to the invention for spectral 
analysis of a wide band signal; 

Figure 2a is a diagram illustrating correlated sub-band signals using the method shown in 
20 Figs, la and lb; 

Figure 2b is a diagram illustrating a comparison of a correlated wide-band spectrum using 
the sub-band correlation method according to the invention and a full band correlation 
method of the prior art; 

25 

Figure 3a is a diagram illustrating percent difference between full band and sub-band 
correlation for identical noise generator seeds; 

Figure 3b is a diagram illustrating percent difference between full band and sub-band 
30 correlation for different noise generator seeds; 

Figure 4a is a diagram illustrating phase difference between full band and sub-band 
correlation for identical noise generator seeds; 
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Figure 4b is a diagram illustrating phase difference between full band and sub-band 
correlation for different noise generator seeds; 

Figure 5 is a is a simplified block diagram of a method according to the invention for 
5 spectral analysis of a wide-band signal; 

Figure 6a is a simplified block diagram of a method according to the invention for 
spectral analysis of wide band signals from a phased array system; and, 

10 Figure 6b is a simplified flow diagram of a method according to the invention for spectral 
analysis of wide band signals from a phased array system. 

Detailed Description of Preferred Embodiments 

15 The wideband interferometric digital architecture WIDAR according to the 

present invention provides efficient wide-band digital signal processing at lower, more 
cost effective system clock rates using parallel processing techniques. Further attributes 
of the invention are high spectral resolution on wide-bands and the ability to quickly 
provide a much higher spectral resolution on smaller arbitrary sub-bands using the same 

20 data. Although, the WIDAR technique is generally applied to signal processing in real- 
time, it is not specifically required for the signals to be real-time. For example, since the 
WIDAR technique reduces the number of operations required it can make non-real time 
analysis on a general purpose computer feasible for some applications. 

25 In the following three applications of the invention a cross-power analyzer, an 

auto-power analyzer and a phased array system will be disclosed. The mathematical 
background of the invention will be explained with respect to the cross-power analyzer. It 
will become evident to a person of skill in the art that the restriction to the cross-power 
analyzer does not unply a loss of generality and the same mathematical background can 

30 be easily applied to other applications. 

Referring to Figs, la and lb a simplified block diagram as well as a simplified 
flow diagram of a method for cross-power analysis according to the present invention is 
shown. A respective system 100 according to the present invention is illustrated in Fig. 

7 
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Ic, showing processing means for processing one analog wide-band input signal. For the 
cross-spectral signal analysis two analog wide-band input signals and Y^^xq 
provided to two systems 100, respectively, each system 100 for processing one signal. 
Final cross-correlation is then performed on one of the two systems. 

5 

In a first processing step the center frequency of the signal is shifted by a 
small fraction of its bandwidth s using an analog mbcer and a local oscillator. The wide- 
band signal is sampled and digitized with a sample rate of at least twice the bandwidth 
of the wide-band signal (Nyquist rate). Optionally, a digital single-sideband mixer 

1 0 may be used to perform the frequency shift digitally after the A/D conversion of the 

signal. The sampled signal is de-multiplexed into iV parallel sample streams for parallel 
processing. Although the step of de-multiplexing is generally used in applications of this 
method it is not strictly required. The parallel sample streams are sent in parallel to 
each of iV digital filters. The bandwidth of each digital filter is approximately MN of the 

1 5 bandwidth of the wide-band signal. Each of the digital filters may have a different tap 
weight. Preferably, digital Finite Impulse Response (FIR) filters such as poly-phase 
decimating FIR filters are used. The sample stream from each digital filter is decimated 
by a factor of , i.e. only every iV* sample is retained and the others are discarded. The 
decimated sample streams are called sub-band signals. In a following step the resolution 

20 of the sub-band signals is reduced by re-scaling and truncating the sub-band signals in 
order to reduce downstream processing load. This step of re-quantization is an optional 
but not necessary feature of the processing. In a following step the sub-band signals are 
phase rotated by phase £■ using a digital phase rotator. Two output signal streams being in 
quadrature with one another are produced for the sample stream from each FIR filter. The 

25 pair of signals is called the de-rotated sub-band signal. TTiere are// de-rotated sub-band 
X signals. 

The processing of the 7 signal is analog, resulting in de-rotated sub-band 
r signals, except that the tap-weights of the FIR filters may be different 

30 

Each of the N de-rotated sub-band A" and F signals are cross-correlated sub- 
band by sub-band. The cross-correlators delay one of the X or the 7 signals with respect 
to the other in a series of delay intervals. The delay intervals are equivalent to one 

8 
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sample-interval of a sub-banci signal. At each of these delays the correlator forms the 
product of the X and 7 signa s in a multiplier. The output of each multiplier is summed 
over a period of time - integration time - producing a sub-band cross-correlation result. 
Each sub-band cross-correlation result is transformed into a cross-spectrum result by 
5 means of a Fourier Transform, which is performed after each integration time. The cross- 
correlators are implemented in system 100 as a software version executed by the N 
processors. Fig. 2a shows an example of a cross-spectrum result for N = \6 sub-bands. 
The sub-band cross-spectra are band flattened and scaled by a cross spectrum weighting 
function. The wide-band spectrum is then obtained by assembling the band flattened and 
10 scaled sub-band cross-spectra side-by side as shown in Fig. 2b. 

Optionally, the sub-band signals may be phase rotated at any time before, during 
or after tiie cross-correlation. 

15 Further optionally, digital sub-sample delay interpolation is used to provide 

precision delay intervals to the sub-band signals. It provides fully digital —: sample delay 

interpolation for delay the delay intervals. After initial quantization, digital delay tracking 
to ± 05 samples of delay is performed with simple integral delays. This imposes a phase 

slope across the full band that is 0 at DC and varies between +^ at the highest 

20 frequency. After filtering and decimation the phase excursion in each sub-band is only 

± samples of delay) plus a phase offset that changes with delay. The phase 

offset in each sub-band can be removed with its fringe stopper. If iV = 16, then 
1 

— sample delay interpolation has been achieved. 
16 

25 Alternatively, the correlator is implemented as a FX correlator. In this correlator, 

the de-rotated sub-band signals are Fourier transformed in real time to the frequency 
domain using, for example. Fast Fourier Transformation (FFT). The real-time spectral 
points are then complex cross-multiplied and tune-averaged to yield the complex cross- 
power spectrum. 

30 
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A digital FIR filter has a predictable amplitude and phase response. The amplitude 
response is the Fourier Transform of its tap coefficients provided the internal word length 
of the realized filter is sufficiently long. For a cosine symmetric FIR filter phase is linear. 
If the integer sample delay through the filter is removed, phase is essentially zero down to 
5 the -15 dB point in the transition band of a bandpass filter. These are the important 
qualities of the FIR filters for the present invention. 

If the sampled signal is de-multiplexed by a factor of N decimation of the output 
of each FIR filter by the same factor (or 2* • N;k = 0,1,...) yields an efficient parallel 
1 0 architecture that has to operate only at the decimated sample rate. Normally, k is set to be 
zero but can be non-zero if it is desired to fiirther split the band into smaller sub-bands 
than produced by the de-muhiplexing. For simplicity, k is set to be zero for the rest of the 

disclosure. If a FIR filter has tap coefficients producing a bandpass of -jj- of the full band 

aligned on sub-band boundaries of ^ ( " = 0,\,...,N - 1 ), decimation by a factor of 

1 5 N yields a baseband signal with bandwidth and sample rate reduced by a factor of , as 
is shown in Crochiere, R. E., Rabiner, T. R. "Multirate Digital Signal Processing", 
Prentice Hall, New Jersey, 1983. If n is odd the signal has an opposite frequency sense 
and is corrected by changing the sign of alternate samples. Any finite-length FIR filter 
has a finite transition band. After decimation, the transition band outside the sub-band 

20 boundaries will suffer aliasing, which causes signals to falsely appear as aliased signals 
within the sub-band. Shifting all spectral features in the wide-band spectrum with a 
frequency shift ^ prevents false correlation of the aliased signals. £-is an arbitrary, but 
small, fi^uency shift., However, shas to be large enough to cause de-correlation in less 
time than the integration time. The aliased signals are quickly suppressed in the output of 

25 the correlator depending upon the length of the integration time and the value of s . The 
rate of suppression is VilsT), wherein T is the integration time. Therefore, false 
correlated signals can be suppressed by an arbitrary amount. If greater suppression is 
needed, the value of s can be varied in a quasi-random way during an integration period. 
If this frequency shift technique is not applied aliased signals will appear and in some 

30 applications aliased "quantization noise" will also appear in the wide-band spectrum. 
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In the following it will be described how correlated sub-band data are corrected to 
allow concatenation with other corrected sub-band data in order to obtain a wide-band 
spectrum. Further details are disclosed by the inventors in Carlson, B.R., Dewdney, B.E. 
"Efficient wideband digital correlation", Electronics Letters, lEE, Vol. 36 No. 1 1, pp987, 
5 May 25, 2000. 

The correction of the correlated sub-band data can be expressed as: 

Pnf = Ponf ^XYFlRr ' "T ^XYi^alr (0 

10 

p„f is a final corrected cross-power spectral point at frequency bin / of the n* sub-band 
from signals X and Y . p^j^ is a cross-power spectral point produced by the sub-band 
correlator (Fourier Transform of correlation of correlation coefficients normalized to «'* 

sub-band X and Y re-quantizer output levels). -— ^ is a sub-band scaling term and 

Pxir 

15 scales each sub-band correlator output to its correct amplitude relative to other sub-bands 
and also to the initial quantizer output. Px^FiRr a gain differential compensation term and 
removes any overall gain differences in sub-band FIR filters inherently present in the 

^XYFJmdeal^ . , , , • , 

denominator of the previous term. — is a bandshape correction term and corrects 

for any sub-band FIR filter shape that is not the same as the ideal filter shape. Finally, 
20 PxTideor ^ wldc-band power gain term that simply scales the output to the desired wide- 
band ideal filter gain. 

^xY - yj{^l){yl) *® geometric mean of the power into the n'" sub-band Xand 
Y re-quantizers. Each power must be measured. The FIR filter tap coefficients across all 
25 sub-bands must be explicitly normalized to a common denominator, usually their 

j_ 

maximum, for (x^) and (3^^) to be meaningful. Pj^^ - [n^Jo Y the geometric mean 
of the power into all re-quantizers. Pxyfir-, o^FXY I is the calculated relative gain 

11 
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of the Xmd F FIR filters compared to the ideal (required) filters, wherein 
^Fxr, - PxYFiR„ I PxfHRideal, • PxYFJE, IS the total cross-power of the sub-band filter 
including aliased power. PxrFm«uai, is the total cross-power of the ideal sub-band filter, 
which has no aliasing. S^^ij^ is the cross-power sub-band filter shape - the product of the 
5 Fourier Transform of the n'" X and 7 sub-band FIR filter tap coefficients. S^^^ . , is 
the desired «'* sub-band filter shape. Pj^y,^^^ is the total power gain of the ideal filter 
fiinction. 

In the following equation (1) will be simplified in order to provide some insight 
10 into the processing of the sub-band correlators. For simplicity, it is useful to re-write 
equation (I) to consider only total power corrections in each sub-band and, for now, 
ignore the individual spectral corrections S^„^^,^ ^ S^fih^ . Equation (1) can thus be re- 
written: 



• • PxYFlRr ■ PxY^r (2) 



In equation (2) p„„ is the normalized sub-band correlation coefficient: 




20 The above equation assumes that there is a linear relationship between the raw correlator 
output {rj) and the normalized correlator output . For coarsely quantized signals, this 
will not generally be true, however the raw correlator output can be normalized and 
linearized using a suitable algorithm. In fact, 4-bit quantization is nearly linear and there 
is little error if the above equation is used directly. 

25 

The second term in equation (2) is the sub-band normalization to the input power 
into the re-quantizers and becomes in its expanded form: 



12 
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where there are N sub-bands. Note that the denominator in the above expression 
is the geometric mean of the power in each sub-band before re-quantization. 

Ahematively, the arithmetic mean [—j; • [S^roV^^X^^)]) "^ed as long as 
PxypfRj. is calculated in a similar fashion. The arithmetic mean generally yields an answer 
that is the same as a full-band correlator, while the geometric mean generally yields a 
more correct answer in the presence of high dynamic range signals. 

The term PxmRr equation (2) is the geometric mean of the calculated cross- 
power gain of all actual sub-band FIR filters relative to the calculated cross-power gain of 
an ideal filter: 



15 



where: 



For an ideal filter with a flat bandpass, the gain is since only of the power 

in the total band is passed by the sub-band filter. The gain Gp„of the real filter is some 
factor Kp„ from the ideal and is: 
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as stated above. 



The final term in equation (2) PxY,deair is the total cross power gain of the ideal 
filter transfer functions. This must be calculated as either the geometric mean or 
arithmetic mean of the power from each sub-band's ideal filter function such that it 
cancels out the gain inherent in the measurement and subsequent calculation of .A 
10 more exact expression for P^,^^^,^ (assuming a geometric mean calculation) is: 



P -(^Ys ] p -flT. 1 

being the number of frequency points in the ideal filter sub-band, and Sxp,R,a,^,^ is the 
ideal filter amplitude at each spectral point / in sub-band « of the X signal (similarly for 
Y). 

Substituting the above expressions into equation (2) and, for now, omitting the 
term yields: 
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wherein G^„„is the re-quantizer power gain of the sub-band « of the Z signal, and 



similarly for the 7 signal. 

10 The square roots of the gains in the above equations are due to the chosen 

convention that the gains are expressed in terms of power and the signals are in terms of 
voltage. Additionally, each decimator has a gain of Gp = N since decimation reduces the 

denominator in the discrete-time calculation of {xl) compared to {x^) in the frequency 
domain convolution of the bandpass signal with the decimation function. Substituting 
15 these expressions into equation (2) yields: 



Since the re-quantizer gains G^^and G^^„as well as the filter gain deviations from ideal 
20 Kpj^„ and Kjry„ are not statistically variant, the above expression becomes: 



which after including Pxr^eair simplifies to: 



(3) 
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From equation (3) it is obvious that the sub-band correlator output is effectively 
normalized to the power out of the initial re-quantizer just like a full band correlator 
would be and that the correct p„ can be obtained if the re-quantizer gains are known. The 
re-quantization gains are dynamic and ensure that near optimum quantization occurs. 
More insight is obtained by considering the response of equation (3) to changes in the 
system. Since the raw sub-band correlator output (jcj?„) is inherently scaled to the power 
into the re-quantizer, any change in that power level due to an increase in uncorrelated 
signal (for example, by using a filter with more aliased power in its particular sub-band) 
will result in a reduction in (;ej)„) . The re-quantizer gain will drop a compensating 
amount and p„ will remain constant. Therefore, the correlated amplitude is independent of 
transition band aliasing. Additionally, if the ideal filter gain increases (for example, if 
more gain is wanted) without adding any more uncorrelated noise, the re-quantizer gains 
will drop and p„ will remain constant. In this case the P^^,^^ term ensures that the final 
output will change with the filter gain change to yield the correct result. 

In the previous section a simplified equation considering the total power in each 
sub-band was derived. Now, the second last term of equation (1) will be included to 
consider its effect. Equation (3) including this term is: 



^"■^ U^^\L-?\' Jg G S W 



wherein the cross spectral power at each spectral point / within each sub-band n is 
considered. 



The band shape correction term: 
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is a result calculated from tht- Fourier Transform of the implemented FIR filter tap 
coefficients and from the required ideal FIR filter frequency response, not including 
aliasing and just using the sp-ictral points / within a particular sub-band. Since it is 
generally not possible for the implemented FIR filter response to be precisely equal to the 
5 ideal FIR filter response, this term provides a correction that can be applied - post 

correlation- to each spectral point to correct the difference- For best interferometer array 
phase closure performance all of the spectral corrections should be close to unity or at 
least symmetric about their respective sub-band center frequencies. 

'0 In the following the effect of the spectral corrections in equation (4) on the 

response of the system will be considered. For example, if S^^jj.^^ is increased by a factor 
of 2 at each frequency bin / (a purposely bad fit of the actual filter response to the ideal 
filter response S^pj^j^,_^ ), the input power to the re-quantizer will increase by a factor of 
2. The re-quantizer gains ( G^„and G,^) will respond by dropping by a factor of 2 and 

'5 A/ will remam constant. This is because simply the level into the re-quantizer is changed. 
The same result will be obtained even if the difference between the actual FIR filter 
amplitude and the ideal FIR filter amplitude is not the same at every frequency bin 
/ »Pon/ always scales to the total power into the re-quantizer, is corrected by a measured 
concomitant amount, and then is corrected at each frequency bin by a calculated ratio. 
20 The ability to correct for virtually any FIR filter fit is an important attribute since it is 
generally true that the FIR filter fit to the required ideal FIR filter shape will be different 
in each frequency bin. 

For comparing the performance of the WIDAR with a fiill band correlator the 
25 ideal filter has unity gain across the sub-band and is zero elsewhere. Figs. 2a and 2b show 
simulation results of a 2 GHz broadband signal containing a strong but relatively 
narrowband signal simulating an astronomical maser. N = \6 was chosen with 4-bit 
initial quantization and 4-bit re-quantization yielding a 250 MHz sub-band sample rate. 
1023-tap FIR filters were used with a sub-band boundary cutoff of -1.2 dB, a flat 
30 passband, and a reject band attenuation of -50 dB. The narrowband signal is slightly 

offset from a sub-band boundary. The bottom trace in Fig. 2b is the full band correlation, 
and the top trace is the sub-band correlation comprising N = \6 concatenated sub-bands, 
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which are shown in Fig. 2a. Three-level fringe stopping was used in the WIDAR sub- 
band correlators and a frequency shift of 24 kHz was introduced in both cases. The full 
band correlation amplitude is reduced because it is inherently normalized to the total 
power into the quantizers, a significant fraction of which is the narrowband signal in this 
5 case. The correlation of the WIDAR sub-band correlator according to the invention was 
set to the correct broadband amplitude by normalizing to one or more sub-band's Pxy„ 
clear of narrowband signals - since only the correlation of a clear band yields the correct 
broadband result. For example, P^^ replaces P^^ and isT^^^ replaces P^^,,^^ in equation 
(1). Fig. 3a shows the results of a channel - by - channel comparison of the WIDAR sub- 

10 band correlator according to the invention and the full band correlator using identical 

noise. As expected, there are peaks at the sub-band boundaries due to aliased independent 
noise from the transition band. These peaks are a very small fraction of the actual noise 
levels. A statistical comparison of noise performance is obtained using independent noise 
as shown in Fig. 3b. Figs. 4a and 4b illustrate the phase comparison for the same two 

15 cases. The sub-band correlation broadband level was found to be within -0.7% of the 
expected level (0.5). This comparison clearly illustrates that the aliased signal 
decorrelates. 

Cross-power analyzers are employed when it is desired to obtain both amplitude 
20 and phase between two signals, which have some common component, i. e. which were 
generated by some common artificial or natural process. In radio astronomy cross power 
analysis along with sufficient averaging time allows very weak signals from distant 
astronomical sources to be extracted from strong independent system noise. Amplitude 
and ph^e from many different antennae are then used to build an image of the radio 
25 source. Other applications for cross-power analyzers include network analysis and 

aperture synthesis radar. In network analysis an artificial signal is generated, transmitted 
through the tested object and is then cross-correlated with the original signal. This yields 
the amplitude and phase response as a function of the frequency of the tested object. The 
tested object could be a circuit or some "black box" responding in different ways at 
30 different frequencies. In aperture synthesis radar an object is irradiated by a high power 
microwave transmitter in order to image the object. One complication with this technique 
is that microwave reflections from many different parts of the object may be coherent 
rather than incoherent making it difficult to reconstruct an image. Here, high spectral 
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resolution on wide-bands would significantly improve the imaging process. In all these 
applications the WIDAR technique according to the invention allows sampler frequencies 
to be much higher than the downstream hardware can process in real time providing a 
significantly higher resolution. 

5 

Referring to Fig. 5 a simplified block diagram of a method for auto-power 
analyzing a signal according to the invention is shown. The method for auto-power 
analysis is a variation of the cross-power analysis method, shown in Figs, la and lb. 
Here, instead of two input signals Xand Yone input signal is split into two. Each of the 
10 two split signals is then sampled independently and cross-correlated as disclosed above 
with respect to the cross-power analyzer. 

Auto-power analyzers are used for spectrum analysis. Digital auto-correlation 
spectrometer allow rapid acquisition of all spectra in a band under consideration but 

15 suffer from limited dynamic range and absolute calibration uncertainty in the presence of 
high dynamic range, narrowband signals. The dynamic range is limited by quantization 
noise and calibration uncertainty arises since the amplitude of the output spectrum is 
always scaled to the power into the quantizer. Therefore, if input signals are changing in 
one part of the band, the entire spectrum amplitude will change. Tliese problems can be 

20 overcome using the WIDAR technique according to the invention. 

Figs. 6a and 6b are diagrams illustrating a method for phased array signal analysis 
according to the invention using the WIDAR technique. A phased array system consists 
of a number of antennas, each with its own receiver. The combined beam of the phased 
array can be steered by electronically changing the delay in each receiver path before 

25 adding the signals. Steering range is determined by the beam of each antenna and steering 
precision is determmed by the accuracy of the delay inserted into each receiver path. 
WIDAR allows phased array applications requiring wide bandwidths, high spectral 
resolution, high dynamic range, and high precision beam steering. Digital sub-sample 
delay interpolation allows insertion of precision delay in receiver paths. It provides fully 

30 digital ^ sample delay interpolation for antenna delay compensation. After initial 

quantization, digital delay tracking to + 0.5 samples of delay is performed with simple 
uitegral delays. 
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The first processing steps of each of the antenna signals are similar to the cross- 
power analyzer shown in Fig. L After re-quantization the sub-band signals are fed to a 
complex mixer followed by 90" phase shift and addition fed with sub-bands provided by 
same digital filters for forming multiple beams. No correlator is shown in Fig. 6. If auto- 
5 correlation is required to allow sub-band spectra to be seamlessly stitched together, then a 
separate mixer, 90* phase shifter, and adder is required to produce the two outputs, which 
have to be provided to the WIDAR auto-correlator. If the phased sub-bands are to be 
cross-correlated with other antennae signals then tibe sub-band outputs would be provided 
to the cross-correlator resulting in cross power spectra with no aliased components. Some 
10 examples of phased array applications are phased array radar, remote sensing phased 
array receivers, and wideband phased antennae. In all these applications the WIDAR 
technique provides wide instantaneous bands and high spectral resolution as well as 
precision beam steering. 

15 Numerous other embodiments of the invention will be apparent to persons skilled 

in the art without departing fi-om the spirit and scope of the invention as defined in the 
appended claims. 
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Claims 

What is claimed is: 

5 LA method for real-time digital spectral analysis of wide-band signals comprising the 
steps of: 

receiving a wide-band signal; 

shifting the center frequency of the wide-band signal by a small fraction e of its 
bandwidth; 

10 sampling and digitizing the wide-band signal; 

processing the digitized wide-band signal using a digital filter; and, 
decimating the digitally filtered wide-band signal. 

2. A method for real-time digital spectral analysis of wide-band signals comprising the 
15 steps of: 

receiving a wide-band signal; 

shifting the center frequency of the wide-band signal by a small fraction fof its 
bandwidth; 

sampling and digitizing the wide-band signal; 

20 de-multiplexing the digitized wide-band signal into N parallel sample streams; 

processing the N parallel sample streams in parallel using AT digital FIR filters; and, 
determining 2'' ■N;k = 0,1,... sub-band signals by decimating the sample stream from 
each digital FIR filter by a factor of 2* • JV;it = 0,1,..., wherein only every 
2*-A^'*;>t = 0,l,... 

25 sample is retained and the others are discarded. 

3. A method for real-time digital spectral analysis of wide-band signals as defined in 
claim 2, wherein the wide-band signal is sampled at a sample rate of at least twice the 
bandwidth of the wide-band signal, 

30 

4. A method for real-time digital spectral analysis of wide-band signals as defined in 
claim 2, wherein each of the N digital FIR filters has a different tap-weight. 
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5. A method for real-time digital spectral analysis of wide-band signals as defined in 
claim 4, wherein each digital FIR filter is a cosine symmetric digital FIR filter having a 
linear phase. 

5 6. A method for real-time digital spectral analysis of wide-band signals as defined in 
claim 5, wherein the bandwidth of each digital FIR filter is approximately 1/ N of the 
bandwidth of the wide-band signal. 

7. A method for real-time digital spectral analysis of wide-band signals as defined in 
10 claim 2, comprising the step of re-quantization by re-scaling and truncating the 

2* • N;k = 0,l,...sub-band signals in order to reduce downstream processing load. 

8. A method for real-time digital spectral analysis of wide-band signals as defined in 
claim 2, comprising the step of phase rotating the 2^ ■N;k=^ 0,1, „. sub-band signals by 

15 phase £ using a digital ph^e rotator producing a de-rotated sub-band signal. 

9. A method for cross-correlating de-rotated sub-band signals sub-band by sub-band, the 
method comprising the steps of: 

receiving 2* - Nik = 0,1,... pairs of first and second de-rotated sub-band signals 
20 at2* • N\k = 0,],... cross-correlators, wherein each pair is received at a different cross- 
correlator of the2* • A'^;^: = 0,1,... cross-correlators; 

delaying one of the first and second de-rotated sub-band signals with respect to the other 
in a series of delay intervals at each of the 2* • JV;^ = 0,1,... cross-correlators; 
forming the product of the first and the second de-rotated sub-band signals at each of the 
25 delay intervals at each of the 2* • JV;>t = 0,1,... cross-correlators; 

producing a sub-band cross-correlation result at each of the 2* •N;k = 0,1,... cross- 
correlators by summing the products over a period of tune; 

transforming each sub-band cross-correlation result at each of the 2* ■ N;k = 0,1,... cross- 
correlators by means of a Fourier Transform into a cross-spectrum result; and, 
30 correcting each cross-power spectral point of each sub-band cross-spectrum result with a 
sub-band scaling term, a gain differential compensation term, a bandshape correction 
term and a wide-band power gain term. 
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10. A method for cross-correlation de-rotated sub-band signals sub-band fay sub-band as 
defined in claim 9, comprising digital sub-sample delay interpolation of the de-rotated 
sub-band signals. 

5 

11. A method for cross-correlating de-rotated sub-band signals sub-band by sub-band as 
defined in claim 9, wherein the delay intervals are equivalent to one sample-interval of a 
sub-band signal. 

10 12. A method for cross-correlating de-rotated sub-band signals sub-band by sub-band as 
defined in claim 9, wherein tiie Fourier Transform is performed after each period of time 
of summing the products. 

13. A method for cross-correlating de-rotated sub-band signals sub-band by sub-band as 
15 defined in claim 9, comprising the step of band flattening each of the 2* • iV; A = 0,1,... 

sub-band cross-spectra. 

14. A method for cross-correlating de-rotated sub-band signals sub-band by sub-band as 
defined in claim 9, comprising the step of scaling each of the 2* • iV;^ = 0,1,.,. sub-band 

20 cross-spectra by a cross spectrum weighting function. 

15. A method for cross-correlating de-rotated sub-band signals sub-band by sub-band as 
defined in claim 9, comprising the step of concatenating the 2* • N;k = 0,1,... sub-band 
cross-spectra to obtain a wide-band spectrum. 

25 

16. A method for cross-correlating de-rotated sub-band signals sub-band by sub-band, the 
method comprising the steps of: 

receiving 2* - N;k = 0,1,,.. pairs of first and second de-rotated sub-band signals 
at2* • iV; A: = 0,1,... cross-correlators, wherein each pair is received at a different cross- 
30 correlator of the 2* • N;k = 0,1,... cross-correlators; 
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transforming each pair of first and second de-rotated sub-band signals at each of 
the2* • N;k = 0,1,... cross-correlators by means of a Fourier Transform into frequency 
domain; 

complex cross-multiplying the Fourier transformed first and second de-rotated sub-band 
5 signals at each of the 2* • 7/; A; = 0,1,. cross-correlators; and, 

time-averaging the cross-multiplied first and second de-rotated sub-band signals. 

17. A method for real-time digital spectral analysis of wide-band signals comprising the 
steps of: 

10 receiving a first and a second wide-band signal; 

shifting the center frequency of each of the first and the second wide-band signal by a 

small fitiction ^jand , respectively, of its bandwidth; 

sampling and digitizing the first and the second wide-band signal; 

de-multiplexing each of the digitized first and second wide-band signals into first// 
15 parallel sample streams and second N parallel sample streams; 

processing each of the first and the second N parallel sample streams in parallel using 

2-iV digital FIR filters; 

determining first and second N sub-band signals by decimating the sample stream from 
each digital FIR filter by a factor of N , wherein only every N"' sample is retained and 
20 the otiiers are discarded; 

re-quantizating the N sub-band signals by re-scaling and truncating in order to reduce 
downstream processing lo^; 

phase rotatuig each of the first and second N sub-band signals by phase and > 
respectively, using a digital phase rotator producing first and second iV^ de-rotated sub- 
25 band signals; 

receiving pairs of the first and second N de-rotated sub-band signals at AT cross- 
correlators, wherein each pair is received at a different cross-correlator of the N cross- 
correlators; 

delaying one of the first and second de-rotated sub-band signals with respect to the other 
30 in a series of delay intervals at each of theiV cross-correlators; 

forming the product of the first and the second de-rotated sub-band signals at each of the 
delay intervals at each of theN cross-correlators; 
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producing a sub-band cross-c«rrelation result at each of theiV cross-correlators by 
summing the products over a period of time; 

transforming each sub-band cross-correlation result at each of theiV^ cross-correlators by 
means of a Fourier Transfona into a cross-spectrum result; 
5 correcting each cross-power spectral point of each sub-band cross-spectrum result with a 
sub-band scaling term, a gain differentia! compensation term, a bandshape correction 
term and a wide-band power gain term; and, 

concatenating the N sub-band cross-spectra to obtain a wide-band spectrum. 

10 1 8, A method for real-time digital spectral analysis of wide-band signals 2S defined in 
claim 17, wherein the frequency shifts s^snd are arbitrary. 

19. A method for real-time digital spectral analysis of wide-band signals as defined in 
claim 1 8, wherein the frequency shifts £, and are varied in a quasi-random way during 

15 a period of time of summing the products. 

20. A method for real-time digital spectral analysis of wide-band signals as defined in 
claim 1 7, wherein the first and the second wide-band signal are a same wide-band signal. 

20 21 . A method for real-time digital spectral analysis of wide-band signals of a phased array 
system comprising the steps of: 

receiving M wide-band signals from the phased array system; 
shifting the center frequency of each of the M wide-band signals by a small fraction 
fii, £2 ,...,£;„, respectively, of its bandwidth; 
25 sampling and digitizing the M wide-band signals; 

de-multiplexing each of tiie digitized M wide-band signals into M-N parallel sample 
streams; 

processing each of the M-'N parallel sample streams in parallel using M-N digital FIR 
filters; 

30 determining M- N sub-band signals by decimating the sample stream from each digital 
FIR filter by a factor of N , wherein only every N''' sample is retained and the others are 
discarded; 

re-quantizing the M-N sub-band signals by re-scaling and truncating; 
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complex mixing each of the M- N re-quantized sub-band signals; 

90° phase shifting one of two components of each of the complex mixed M- N sub-band 

signals; and, 

forming multiple beams by adding same sub-band components of the M- N sub-bands, 
5 the same sub-band components being provided by same digital FIR filters. 

22. A method for real-time digital spectral analysis of wide-band signals as defined in 
claim 21, comprising digital sub-sample delay interpolation of the digitized M wide- 
band signals. 

10 

23. A system for real-time digital spectral analysis of wide-band signals comprising: 
a port for receiving a wide-band signal; 

a frequency shifter for shifting the center fi-equency of the wide-band signal by a small 
fraction £r of its bandwidth; 
15 an A/D converter for sampling and digitizing the wide-band signal; 

a de-multiplexer for de-multiplexing the digitized wide-band signal into N parallel 
sample streams; and, 

N processors, each processor for processing one of the N parallel sample streams by 
digitally FIR filtering and decimating the sample stream in order to determine a sub-band 



24. A system for real-time digital spectral analysis of wide-band signals as defined in 
claim 23, wherein the frequency shifter the center frequency comprise an analog mixer 
and a local oscillator. 



25. A system for real-time digital spectral analysis of wide-band signals as defined in 
claim 23, wherein the fi-equency shifter the center frequency comprise a digital single- 
sideband mixer. 

30 26. A system for real-time digital spectral analysis of wide-band signals as defined in 
claim 23, wherein each of the N processors re-quantizes the sub-band signal. 



20 



signal. 
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27. A system for real-time digital spectral analysis of wide-band signals as defmed in 
claim 23, wherein each of the N processors digitally phase rotates the sub-band signal by 
phase-rate e, 

5 28. A system for real-time digital spectral analysis of wide-band signals as defined in 
claim 27, wherein each of the N processors digitally cross-correlates the sub-band signal 
with a respective second sub-band signal. 



10 
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receiving a first and a second wide-band signal 








shifting the center frequency of each of the first and the second 

wide-band signal by a small fraction and ^2 , respectively, of 
its bandwidth 






r 




sampling and digitizing the first and the second wide-band signal 



I 



de-multiplexing each of the digitized first and second wide-band 
signals into first N parallel sample streams and second at parallel 
sample streams 

processing each of the first and the second N parallel sample 
streams in parallel using 2-n digital FIR filters 

i ~" 

determining first and second sub-band signals by decimating the 
sample stream from each FIR filter by a factor of N , wherein only 
every N'^ sample is retained and the others are d iscarded 

re-quantizing the N sub-band signals by re-scaling and 
truncatmg in order to reduce downstream processing load 

i 

TO FIG IB (cont.) 
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FROM FIG. IB 

n 


phase rotating each of the first and second n sub-band signals 
by phase-rate and ^2 , respectively, using a digital phase 
rotator producing first and second de-rotated sub-band signals 






receiving pairs of the first and second n de-rotated sub-band 
signals at cross-correlators, wherein each pair is received at a 
different cross-correlator of thcA^ cross-correlators 




X 


delaying one of the first and second de-rotated sub-band signals 
with respect to the other m a series of delay intervals at each of 
theiv cross-correlators 





forming the product of the first and the second de-rotated sub- 
band signals at each of the delay intervals at each of the N cross- 
correlators ^ 

i 

producmg a sub-band cross-correlation result at each of the N 
cross-correlators by summing the products over a period of time 

i 

transforming each sub-band cross-correlation result at each of 
theiv cross-correlators by means of a Fourier Transform into a 
cross-spectrum result 



assembling theiv sub-band cross-spectra side by side to obtain a 
wide-band spectrum 

Fig. lb (cont.) 
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Processor for: 
digitally FIR filtering 

decimating for determining sub-band signals 
re-quantizing the sub-band signals 
digitally phase rotating the sub-band signals 
digitally cross-correlating the sub-band signal with a 
respective second sub-band signal 
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receiving M wide-band signals from the phased array system 

i 

shifting the center frequency of each of the M wide-band signals 
by a small fraction 5 ^2 » • • • » , respectively, of its bandwidth 

i ~" 

sampling and digitizing the m wide-band signals 

i 

de-multiplexing each of the digitized M wide-band signals into 
M-N parallel sample streams 

processing each of the m-n parallel sample streams in parallel 
using M-N digital FIR filters 

determining m- n sub-band signals by decimating the sample 
stream from each FIR filter by a factor of n, wherein only every 

N*^ sample is retained and the others are discarded 

re-quantizing the m-n sub-band signals by re-scaling and truncating 

4 

complex mixing each of the m-n re-quantized sub-band signals 
and off-setting their phase to effect digital sub-sample delay 
interpolation 

90° phase shifting one of two components of each of the complex 
mixed m-n sub-band signals 

i 

forming multiple beams by adding same sub-band components of 
the m-n sub-bands, the same sub-band components being provided 
by same digital FIR filters 

Fig. 6b 
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